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Abstract 

We study the heavy-quark contributions to the proton structure function F2(x,Q 2 ) at 
next-to-leading order using compact formulas at small values of Bjorken's x variable. The 
formulas provide a good agreement with the modern HERA data for F^x^Q 2 ). 
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1 Introduction 

In the last year the HI [TJ [2] and ZEUS [31 0] Collaborations at HERA are presented a new 
data on the charm structure function (SF) F£ 0- Moreover, recently HI and ZEUS have been 
demonstrated the preliminary combine data of F£(x, Q 2 ) [5]. 

In the framework of Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) dynamics [6j 
there are two basic methods to study heavy-flavor physics. One of them [7] is based on the 
massless evolution of parton distribution functions (PDFs) and the other one is based on 
the photon-gluon fusion (PGF) process [8]. There are also some interpolating schemes (see 
Ref. [9] and references cited therein). 

In this short paper, we apply compact low-x approximations for the SF F%{x, Q 2 ) obtained 
|10| in the PGF framework at the first two orders of perturbation theory to these new HERA 
experimental data [lj [2j [31 HI [5] . We show a good agreement between experimental data and 
the approach which was found without any additional free parameters. All PDF parameters 
have been fitted earlier from F2(x,Q 2 ) HERA experimental data only. 



2 Approach 

We now present our basic formula for F^ix^Q 2 ) appropriate for small values of x, where 
only the gluon and quark-singlet contributions matter, while the non-singlet contributions 
are negligibly small □ 

F^(x,Q 2 )= J2 E C£ a (x,Q 2 ,^ 2 )0x/:(x,^)« J2 C r 2 \ g {x,Q 2 ^ 2 )® x f g {x^ 2 ) (1) 

a=g,l,l r =+- r =+>- 

1 Thc papers [l][2l|3l|4] contain also the references on the previous data on deep-inelastic (DIS) structure function 
(SF) F| at small x values. 

2 Here and in the following, we suppress the variables m c in the argument lists of the structure and coefficient 
functions for the ease of notation. 



where I generically denotes the light-quark flavors, r = ± labels the usual + and — linear 
combinations of the gluon and quark-singlet contributions, C 7 2 a {x) are the DIS coefficient 
functions, which can be calculated perturbatively in the parton model of QCD, \x is the 
factorization scale and same time the renormalization one appearing in the strong-coupling 
constant a s (fj,), and the symbol <8> denotes convolution according to the usual prescription, 
f(x) ® g(x) = f^(dy/y)f(y)g(x/y). Massive kinematics requires that C\ a = for x > b = 
1/(1 + 4c), where c = m 2 /Q 2 . We take m c to be the solution of m c (m c ) = m c , where m c (/i) is 
defined in the modified minimal-subtraction (MS) scheme. The simplification in the r.h.s. of 
(PQ) is obtained by neglecting the contributions due to incoming light quarks and antiquarks, 
which is justified because they vanish at LO and are numerically suppressed at NLO for 
small values of x. 

Exploiting the low-x asymptotic behavior of f^(x,Q 2 ) [12\ 113] . 

fa( x ,^ 2 ) X -$° 7^s;fa( x ^ 2 ) (hereafter a = q,g), (2) 

where the rise of (x, /U 2 ) as x — > is less than any power of x, Eq. (prj can be rewritten as 
FZ(x,Q 2 )* £ M£ 9 (l + 5 r ,QV) */;(*, /, 2 ), (3) 

r=+,- 

where 

MJ» = j\xx n - 2 Cl g (x) (4) 

is the Mellin transform, which is to be analytically continued from integer values n to real 
values 1 + 8 r . Strictly speaking, the equation ([3]) is correct for non-singular Mellin moments 
M% g (n) at n — > 1. The generalization of ([3]) for the moments containing singularity will be 
done in subsection 4.2. 



3 Gluon density 

As demonstrated in Refs. [141 HERA data for i^O^, Q 2 ) support the modified Bessel-like 
behavior of PDFs at small x values predicted in the framework of the so-called generalized 
double-asymptotic scaling regime. El In this approach, DGLAP dynamics [6] starting at some 
initial value /jLq with flat x distributions: 

xfa(x,/4) = A a, (5) 

where A a are unknown constants to be determined from the data. 

In the framework of the generalized double-asymptotic scaling regime, the main ingredi- 
ents of the results for gluon density include the following0 It is presented in terms of two 
components ("+" and "— ") 

f g (x,v 2 ) = f+(x,fi 2 ) + f-(x,[i 2 ), (6) 

which are obtained from the analytic /^-dependent expressions of the corresponding ("+" 
and "— ") PDF moments. Here, ej is the fractional electric charge of heavy quark i, e = 
(J2i=i e i)// i s the average charge square and / is the number of active quark flavors. 

The small-x asymptotic results for the PDF fg are at the next-to-leading order (NLO) 

m ' 

xf+{x^ 2 ) = Apoi^e-^-^ + Oip), 

xf-(x,p 2 ) = A-e- d - s - D -P + 0{x), (7) 

3 It is a generalization of earlier studies [T5] . 

4 The results for quark densities may be found in Refs. [Til 111) . 



where I u are the modified Bessel functions, 



D± = d±± - —d± 

Po 



(8) 



and similar for D + and D + , 
At = 
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with £3 and £2 are Euler functions. /3q and f}\ are first two coefficients of QCD /3-function 
The variables s, p, a, and p are given by 

a s {po 



In 



a s (p) 



p = a s (p Q ) - a s (p), a 



s + D + p \ lux, p 



a 



21n(l/x) 
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The transformation to the LO is simple: in Eqs. (|7|), ([9]) and (llip we should put 
p = a s (Q) = a s (Qa) = and replace the variable s by its LO approximation: s LO = 
Ma™(p )/a^°(p)). 



4 Wilson coefficients 



One has (1) = (1) if M^giri) are devoid of singularities in the limit 5 r —> 0, as we 
assume for the time being. Such singularities actually occur at NLO, leading to modifications 
to be discussed below. Defining M2 )9 (l) = M^ g (l) and using f g (x,Q 2 ) = J2r=± fg ( x ; Q 2 ) > 
Eq. (J3]) may be simplified to become 

F 2 c (x, Q 2 ) « M 2jfl (l, Q 2 ,p 2 )xf g (x, p 2 ). (12) 

Through NLO, M 2 , g (l,Q 2 , p 2 ) exhibits the structure 

M 2 , 9 (l,Q 2 ,/i 2 ) = e 2 a s (/x) { M^(l, c) + a s (/x) 



Mg(l,c) + Mg(l, C )ln^ 



+ 0{a\ 



(13) 



4.1 LO results 

The LO coefficient function of PGF can be obtained from the QED case [T7f by adjusting 
coupling constants and color factors, and they read [El [19] 



Ci 0) g(x,c) = -2s{[l-4x(2-c)(l-x)]/3- [1 - 2x(l - 2c) + 2x 2 (l - 6c - 4c 2 )]L(/3)}, (14) 



where 



(3(x) 



1 



4cx 



L(/3) = In 



1 + /3 



1-x' v ~ y ""1-/3 
Performing the Mellin transformation [10] in Eq. ([5]), we find 



m2(1,c) = + 2(1 -c)J(c)] 



with 



J(c) = -Vb hit, i 



1-V6 
1 + Vfe' 



1 

1 + 4c' 



(15) 

(16) 
(17) 



4.2 NLO results 



The NLO coefficient functions of PGF are rather lengthy and not published in print; they 
are only available as computer codes [20J. For the purpose of this letter, it is sufficient to 
work in the high-energy regime, defined by x <C 1, where they assume the compact form |21j 



C^ g (x,c) = PR^ g (l,c), (18) 

with 

4>,c) = ^C A [5 + (13-10c)J( C )+6(l- C )/( C )], 42(1,0) = -4C A M 2 ( 2(1 >C ), (19) 
where Ca = N for the color gauge group SU(N), J(c) is defined by Eq. (fTT|) . and 



1(c) = -Vb 



C(2) + - In 2 t - ln(6c) In t + 2 Li 2 (— *) 



(20) 



where t is given in ()1T|) and Li2(xc) = — f}(dy/y) ln(l — ccy) is the dilogarithmic function. 

As already mentioned above (see the end of Section 2), the Mellin transforms of C^ g {x, c) 
exhibit singularities in the limit 5 r — > 0, which lead to modifications in Eqs. ([3j) and (112ft . 
As was shown in Refs. |X3|, [T41 [TT] . the terms involving l/<5 r correspond to singularities of 
the Mellin moments M^gin) at n — ¥ 1 and depend on the exact form of the subasymptotic 
low-x behavior encoded in f g (x,fj?). The modification is simple: 

11 1 1 f 1 dy ~ r 



hjr/ %(^ 2 ), (21) 



6r S r 8 r f' g (x^ 2 ) h V 

where x = x/b. In the generalized double-asymptotic scaling regime, the + and — compo- 
nents of the gluon PDF exhibit the low-x behavior (|7|) . We thus have |10t [TT] 

<5 + />(£) 7 (cr(a;)) J_ a; 



where cr and p are given in (jlip . 

Because the ratio f g (x, Q 2 )/f+(x, Q 2 ) is rather small at the Q 2 values considered, Eq. (fT2j) 
is modified to become 

F 2 c (x, Q 2 ) « M 2)fl (l, fi 2 ,c)xf 9 (x, /i 2 ), (23) 

where M2,g(l, A* 2 ) is obtained from M2 jg (n, p 2 ) by taking the limit n — )• 1 and replacing 
l/(n — 1) — »• l/&f. Consequently, one needs to substitute 

M 2 <2(l,c)->Mg(l,c) (J = 1,2) (24) 
in the NLO part of Eq. ([13]) . Using the identity 

f -/5 (-) W(V)) « i - In(fec) - (25) 

find the Mellin transform © of Eq. {[TED to be El 

i?g(l,c) (j = l,2), (26) 



wc 



Mg(l,c) 



57 - ln(6c) - — 



with R2 g (l,a) (j = 1,2) are given in (|26p . The rise of the NLO terms as x — > is in 
agreement with earlier investigations [24J. 

5 Note, that 6+ determines the behavior of the slope of gluon density (see and also mostly the slope of SF 
Fi . The form ([2"2"j) of 5+ is in good agreement 22 with the corresponding HERA experimental data. 



5 Results 



We are now in a position to explore the phenomenological implications of our results. As 
for our input parameters, we choose m c = 1.25 GeV in agreement with Particle Data Group 
[23j . While the LO result Eq. (|16H is independent of the unphysical mass scale fj,, the NLO 
formula dT3j) does depend on it, due to an incomplete compensation of the /x dependence 
of a(fJL) by the terms proportional to ln(/i 2 /Q 2 ), the residual [i dependence being formally 
beyond NLO. In order to fix the theoretical uncertainty resulting from this, we put fj? = 
Q 2 + 4m 2 , which is the standard scale in heavy quark production. 

The PDF parameters //{j, A q and A g shown in ([5]), (|7|) and Q have been fixed in the fits 
of Fi experimental data. Their values depend on conditions chosen in the fits: the order of 
perturbation theory and the number / of active quarks. 

Below 6-quark threshold, the scheme with / = 4 has been used [TT] in the fits of F<i 
data. Note, that the F% structure function contains F£ as a part. In the fits, the NLO gluon 
density and the LO and NLO quark ones contribute to F^ , as the part of Fi- Then, now in 
PGF scattering the LO coefficient function (I14p corresponds in m — > limit to the standard 
NLO Wilson coefficient (together with the product of the LO anomalous dimension ^ qg and 
ln(m 2 /Q 2 )). It is a general situation, i.e. the coefficient function of PGF scattering at some 
order of perturbation theory corresponds to the standard DIS Wilson coefficient with the one 
step higher order. The reason is following: the standard DIS analysis starts with handbag 
diagram of photon-quark scattering and photon-gluon interaction begins at one-loop level. 

To analyze F| at the LO of PGF process one should take xf a (x, Q 2 ) from the fit of F2 at 
NLO with / = 4. In practice, we use the following parameters, Qq = 0.523 GeV 2 , A g = 0.060 
and A q = 0.844. 

Correspondingly, to analyze PGF process at NLO one needs to know the gluon density 
extracted from F2 data at NNLO, which is not yet known |f| in generalized double-asymptotic 
scaling regime. However, as we can see from the modern global fits [26], the difference 
between NLO and NNLO gluon densities is not so large. So, we can safely apply the NLO 
form (J7D of xf g (x,Q 2 ) to our NLO PGF analysis. 

The results for F% are presented in Fig.l. We can see a good agreement between our 
compact formulas (fl~2l) . (fTBj) and ([26]) and the modern experimental data [TJ [2J [3j [4J, [5] for 
F^x^Q 2 ) structure function. To keep place on Fig.l, we show only the HI [H [2] and the 
combine preliminary Hl-ZEUS [5] data. 

The good agreement between generalized double-asymptotic scaling approach and F2, F£ 
data demonstrates an equal importance of the both parton densities (gluon and sea quark) 
at low x. It is due to the fact that F% relates mostly to the sea quark distribution, while the 
F£ relates mostly to the gluon one. Dropping sea quarks in analyze leads to the different 
gluon densities extracted from F2 or from F% (see, for example, [27] ) . 

6 Conclusions 

We presented a compact formulas for the heavy-flavor contributions to the proton structure 
functions F2 valid through NLO at small values of Bjorken's x variable. Our results agree 
with modern experimental data [U [21 [3j [H [5] well within errors without a free additional 
parameters. In the Q 2 range probed by the HERA data, our NLO predictions agree quite well 
with the LO ones. Since we worked in the fixed-flavor-number scheme, our results are bound 
to break down for Q 2 S> 4m 2 , which manifests itself by appreciable QCD correction factors 
and scale dependence. As is well known, this problem is conveniently solved by adopting the 
variable-flavor-number scheme, which we leave for future work. 

6 The difficulty to extend the analysis [THE] to NNLO is related to the appearance of the pole ~ l/(n— l) 2 in 
the three- loop corrections to the anomalous dimension j gg (see [25] )• The pole ~ l/(n— l) 2 violates the Bessel-like 
solution ([7]) of DGLAP equation for PDFs at low x values with the flat initial condition ([5j. 
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Figure 1: (colored online) The charm structure function F£(x, Q 2 ) evaluated as function of x with 
the LO matrix elements from Eq. (fl~6l) (dashed lines) and with the NLO ones from Eq. fl26|) and 
with the factorization/renormalization scale fi 2 = Q 2 -\-Am 2 (solid lines). The black points and red 
squares correspond to the combine preliminary Hl-ZEUS data [S] and HI data [HE], respectively. 

Acknowledgements 

A.V.K. was supported in part by the Alexander von Humboldt Foundation. The work was 
supported in part by RFBR grant No.l0-02-01259-a. 

Calculations were partially performed on the HPC facility of SISSA/Democritos in Trieste 
and partially on the HPC facility "WIGLAF" of the Department of Physics, University of 
Trento. 

References 

[1] F. D. Aaron et al. [HI Collaboration], Phys. Lett. B 686 (2010) 91 [arXiv:0911.3989l 
[hep-ex]]. 

[2] F. D. Aaron et al. [HI Collaboration], Eur. Phys. J. C 65 (2010) 89 [arXiv:0907.2643l 
[hep-ex]]. 

[3] H. Abramowicz et al. [ZEUS collaboration], arXiv: 1005.3396 [hep-ex]. 

[4] S. Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. C 65 (2010) 65 [arXiv:0904.3487l 
[hep-ex]]. 



[5] K. Lipka [HI Collaboration and ZEUS Collaboration], Nucl. Phys. Proc. Suppl. 191 
(2009) 163; A. M. Cooper-Sarkar. larXiv: 1006.4471 1 [hep-ph]. 

[6] V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438, 675 [Yad. Fiz. 15 
(1972) 781, 1218]; L. N. Lipatov, Sov. J. Nucl. Phys. 20 (1975) 94 [Yad. Fiz. 20 (1974) 
181]; G. Altarelli and G. Parisi, Nucl. Phys. B 126 (1977) 298; Yu. L. Dokshitzer, Sov. 
Phys. JETP 46 (1977) 641 [Zh. Eksp. Teor. Fiz. 73 (1977) 1216]. 

[7] B. A. Kniehl, G. Kramer and M. Spira, Z. Phys. C 76 (1997) 689 
arXiv:hep-ph/9610267|; J. Binnewies, B. A. Kniehl and G. Kramer, Z. Phys. 
C 76 (1997) 677 [arXiv:hep-ph/9702408| ; Phys. Rev. D 58 (1998) 14014 
[arXiv:hep-ph/9712482| ; Phys. Rev. D 58 (1998) 034016 [arXiv:hep-ph/98022311 . 

[8] S. Fri xione, M. L. Man gano, P. Nason and G. Ridolfi, Phys. Lett. B 348 (1995) 633 
[arXiv:hep-ph/9412348]; S. Frixione, P. Nason and G. Ridolfi, Nucl. Phys. B 454 (1995) 
3 [arXiv:hep-ph/9506226] . 

[9] F. I. Olness and W. K. Tung, Nucl. Phys. B 308 (1988) 813; M. A. G. Aivazis , 
F. I. Olness and W. K. Tung, Phys. Rev. D 50 (1994) 3085 |arXiv:hep-ph/9312318| ; 
M. A. G. Aivazis, J. C. Collins, F. I. Olness and W. K. Tung, Phys. Rev. D 50 (1994) 
3102 |arXiv:hep-ph/9312319] ; R. S. Thorne and R. G. Roberts Phys. Rev. D 5 7 (1998) 
6871 [arXiv:hep-ph/9709442| ; P hys. Lett. B 421 (1998 ) 303 [arXiv:hep-ph/9711223| ; 
Eur. Phys. J. C 19 (2001) 339 [arXiv:hep-ph/0010344]; W. K . Tung, S. Kretzer and 
C. Schmidt, J. Phys. G 28 (2002) 983 [arXiv:hep-ph/011024 7|. 

[10] A. Y. Illarionov, B. A. Kniehl and A. V. Kotikov, Phys. Lett. B 663 (2008) 66 
[arXiv:0801. 15021 [hep-ph]]. 

[11] A. Y. Illarionov, A. V. Kotikov and G. Parente Bermudez, Phys. Part. Nucl. 39 (2008) 
307 ]arXiv:hep-ph/0402173| . 

[12] C. Lopez and F. J. Yndurain, Nucl. Phys. B 171 (1980) 231; Nucl. Phys. B 183 (1981) 
157. 

[13] A. V. Kotikov, Phys. Rev. D 49 (1994) 5746; Phys. Atom. Nucl. 59 (1996) 2137 [Yad. 
Fiz. 59 (1996) 2219]. 

[14] A. V. Kotikov and G. Parente, Nucl. Phys. B 549 (1999) 242 [arXiv:hep-ph/9807249] . 

[15] A. De Rujula, S. L. Glashow, H. D. Politzer, S. B. Treiman, F. Wilczek and A. Zee, 
Ph ys. Rev. D 10 (1974) 1649; R. D. Ball and S. Forte, Phys. Lett. B 335 (1994) 77 
[arXiv:hep-ph/9405320 1 ; L. Mankiewicz, A. Saalfeld and T. Weigl, Phys. Lett. B 393 



(1997) 175 arXiv:hep-ph/9612297|. 



[16] A. V. Kotikov, Phys. Part. Nucl. 38 (2007) 1 [Erratum-ibid. 38 (2007) 

[17] V. N. Baier, V. S. Fadin and V. A. Khoze, Sov. Phys. JETP 23 (1966) 104 [Zh. Eksp. 
Teor. Fiz. 50 (1966) 156]; V. G. Zima, Yad. Fiz. 16 (1972) 1051; V. M. Budnev, 
I. F. Ginzburg, G. V. Meledin and V. G. Serbo, Phys. Rept. 15 (1974) 181. 

[18] E. Witten, Nucl. Phys. B 104 (1976) 445; J. P. Leveille and T. J. Weiler, Nucl. Phys. 
B 147 (1979) 147; V. A. Novikov, M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, 
Nucl. Phys. B 136 (1978) 125 [Yad. Fiz. 27 (1978) 771]. 

[19] A. V. Kotikov, A. V. Lipa tov, G. Parente and N. P. Zotov, Eur. Phys. J. C 26 (2002) 
51 | arXiv:hep-ph/0107135] . 

[20] E. Laenen, S. Riemersma, J. Smith and W. L. van Neerven, Nucl. Phys. B 392 (1993) 
162, 229. 

[21] S. Catani, M. Ciafaloni and F. Hautmann, Preprint CERN-Th.6398/92, in Proceeding 
of the Workshop on Physics at HERA (Hamburg, 1991), Vol. 2., p. 690; S. Catani and 
F. Hautmann, Nucl. Phys. B 427 (1994) 475 [arXiv:hep-ph/9405388|; S. Riemersma , 
J. Smith and W. L. van Neerven, Phys. Lett. B 347 (1995) 143 [arXiv:hep-ph/9411431| . 



[22] A. V. Kotikov and G. Parente, J. Exp. Theor. Phys. 97 (2003) 859 [Zh. Eksp. Teor. Fiz. 
97 (2003) 963] [arXiv:hep-ph/0207276] ; G. Cvetic, A. Y. Illarionov, B. A. Kniehl and 
A. V. Kotikov, Phys. Lett. B 679 (2009) 350 [arXiv: 0906 .19251 [hep-ph]]. 

[23] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667 (2008) 1. 

[24] P. Nason, S. Dawson and R. K. Ellis, Nucl. Phys. B 303 (1988) 607. 

[25] V. S. Fadin and L. N. Lipatov, Phys. Lett. B 429 (1998) 127 |arXiv:hep-ph/9802290| ; 
A. V. Kotikov and L. N. Lipatov, Nucl. Phys. B 582 (2000) 19 [arXiv:hep-ph/0004008] ; 
A. Vogt, S. Moch and J. A. M. Vermaseren, Nucl. Phys. B 691 (2004) 129 
[arXiv:hep-ph /0404m] . 

[26] M. Dittmar et al, |arXiv:hep-ph/0511lT9| 

[27] H. Jung, A. V. Kotik ov, A. V. Lipatov and N. P. Zotov, larXiv:0706.3793 [hep-ph]; 
|arXiv:hep-ph/06110"93| 



